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1.  Introduction 

In  this  paper,  a  continuum  model  for  materials  that  can  undergo  martensitic  phase 
transformations  is  developed  and  applied  to  the  study  of  several  problems  that  involve 
such  transformations.  Among  other  things,  the  continuum  model  that  is  developed 
provides  the  correct  material  symmetry  group  for  each  phase  of  the  material,  and 
results  in  the  corresponding  boundary  value  problem  being  in  a  form  that  permit 
direct  linearization,  while  still  retaining  finite  shape  deformations  for  the  martensite 
phases.  All  of  the  problems  that  this  continuum  model  is  applied  to  in  this  paper 
deal  with  the  issue  of  which  phase  or  which  variant  of  martensite  is  preferred  during 
the  growth  process  when  a  boundary  traction  is  applied.'  Among  these  problems  are 
the  case  of  a  uniaxial  tensile  traction  applied  to  a  cylindrical  body,  and  the  case  of 
a  hydrostatic  pressure  applied  to  a  material  that  has  a  finite  sh^  deformation  with 
an  infinitesimal  dilatation. 

The  term  martensitic  phase  transformation  was  originally  given  to  the  diffu* 
sionless  phase  transformation  that  occurs  when  the  high  temperature  austenite  phase 
of  a  steel  is  quenched,  and  the  term  martensite  was  originally  given  to  the  phase 
that  is  created  from  this  solid-solid  phase  transformation.  The  term  martensitic  phase 
transformation  has  since  been  given  to  almost  all  solid-solid  phase  transformations  that 
proceed  by  a  diffusionless  cooperative  movement  of  atoms  at  the  phase  boundary, 
involve  a  change  in  crystal  structure,  have  a  distinct  orientation  relation  between 
the  crystal  lattices  of  the  parent  phase  and  the  product  phase,  have  a  deformation 
associated  with  the  product  phase,  and  have  continuity  of  displacements,  with  possible 
discontinuity  of  strain,  at  the  phase  boundary.^  Similarly,  the  phase  that  is  created 


'  This  continuuni  model  was  also  applied  to  the  study  of  the  terapenture  at  the  inteifaoe  and  the  quasi-static 
motions  of  a  two-phase  thennoelastic  bar  [9],  [10].  and  to  the  study  of  the  longitudinal  fiee  vibrations  of  a  finite 
fixed-free  two-phase  elastic  bar  [9],  [11]. 

*  See  [8]  and  [15]  for  more  detailed  and  comprehensive  disucussions  about  martensitic  phase  itinsfoimatioas. 
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by  a  transformation  that  is  considered  to  be  martensitic  is  sometimes  referred  to 
as  the  martensite  phase  of  the  given  material.  In  this  paper,  when  a  martensitic 
phase  transformation  between  a  high  temperature  phase  and  a  low  temperature 
phase  of  a  material  is  being  considered,  the  high  temperature  phase  will  usually 
be  referred  to  as  the  austenite  phase,  and  the  low  temperatiue  phase  will  usually  be 
referred  to  as  the  martensite  phase,  regardless  of  which  phase  is  being  created  by  the 
martensitic  phase  transformation.  Additionally,  the  martensitic  phase  transformation 
that  creates  the  austenite  phase  may  son^times  be  referred  to  as  the  austenite  phase 
transformation,  and  the  martensitic  phase  transformation  that  creates  the  martensite 
phase  may  sometimes  be  referred  to  as  the  martensite  phase  trarrsformation. 

As  mentioned  above,  the  martensite  phase  is  characterized  by  having  a 
deformation  relative  to  the  tmdeformed  parent  phase.  This  deformation  consists 
mostly  of  the  deformation  that  occurs  solely  from  the  mechanisms  of  the  martensitic 
phase  transformation  that  occur  at  the  phase  boundary.  This  portion  of  the  total 
deformation  corresponds  to  an  unstressed  state  of  the  martensite  and  will  henceforth 
be  referred  to  as  the  shape  deformation  of  the  martensite.^  Additionally,  the 
strain  corresponding  to  the  shape  deformation  will  sometimes  be  referred  to  as 
the  transformation  strain.  The  shape  deformation  is  in  general  a  finite  deformation, 
and  consists  primarily  of  the  deformation  that  would  be  necessary  to  deform  the 
austenite  crystal  lattice  into  the  martensite  crystal  lattice  and,  for  some  materials, 
the  deformation  that  is  necessary  to  maintain  continuity  of  displacements  at  the 
phase  boundary.  The  remaining  portion  of  the  total  deformation  corresponding  to  the 
martensite  is  due  to  the  surrounding  matrix  material  constraining  the  formation  of  the 
shape  deformation  and/or  by  any  applied  boundary  tractions.  These  deformations  are 

^  Becauw  the  shape  defoniiation  conesponds  to  an  unstressed  state  of  the  martensite,  it  can  be  considered  to 
represent  the  undefonned  martensite. 
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usually  infinitesimal  deformations,  even  though  they  may  be  greater  than  the  yield 
strain  of  the  martensite,  and  in  this  paper  they  will  be  considered  to  be  superimposed 
upon  the  finite  shape  deformation.  If  the  martensite  is  stressed,  it  is  due  solely  to 
these  superimposed  deformations. 

The  change  in  crystal  structure  that  occurs  during  a  martensitic  phase  transforma¬ 
tion  may  involve  a  change  in  crystal  symmetry,  as  is  the  case  with  austenite-martensite 
phase  transformations  in  most  materials  that  can  undergo  such  transformations,  or 
it  may  result  in  the  product  phase  having  the  same  type  of  crystal  lattice  as  the 
parent  phase,  but  with  a  difierent  orientation,  as  is  the  case  when  one  variant  of 
martensite  is  transformed  into  another  variant  of  the  same  martensite.  The  symmetry 
of  the  crystal  lattice  of  a  material  is  reflected  in  the  material  symmetry  group  of  the 
material.  The  material  symmetry  group  of  a  material  restricts  the  functional  form 
of  the  constitutive  equations  of  the  material  (see  [13]).  Thus,  when  a  particle  of 
material  in  its  austenite  phase  is  transformed  to  its  martensite  phase,  the  material 
symmetry  group  of  that  particle  of  material  should  change  accordingly,  and  this 
should  be  reflected  in  the  constitutive  equations  of  that  particle  of  material  both 
before  and  after  the  transformation.  This  change  in  material  symmetry  group  should 
not  only  represent  the  change  in  the  type  of  crystal  lattice,  it  should  also  represent  the 
orientation  relation  between  the  crystal  lattices  of  the  austenite  and  the  martensite. 

In  the  next  few  sections,  the  continuum  model  for  materials  that  can  undergo 
martensitic  phase  transformations  is  developed  and  the  corresponding  field  equations 
and  jump  conditions  are  derived  for  a  purely  mechanical  process.  A  similar  continuum 
model  can  be  developed  and  the  corresponding  field  equations  and  jump  conditions 
can  be  derived  for  a  process  that  is  thermo-mechanical  in  a  similar  manner  (see  [9]). 
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2.  The  Eulerian  Global  Form  of  the  Balance  Laws 

It  is  assumed  that  the  process  under  consideration  takes  place  at  a  constant, 
uniform  temperature  and  with  no  heat  conduction.  Such  a  process  is  an  isothermal 
and  adiabatic  process  and  is  also  known  as  a  purely  mechanical  process.  It  is 
assumed  that  this  process  takes  place  in  a  time  interval  P  =  The  body  B 

that  is  considered  is  assumed  to  occupy  a  regular  region  R«  at  time  t  €  P  and  it  is 
assumed  that  Rt  is  a  subset  of  the  three-dimensional  Euclidean  space  E3.  A  point 
or  the  position  vector  of  a  point  in  Rt  is  denoted  by  y.  The  traction  on  the  surface 
with  unit  normal  n(y,t)  is  denoted  by  t(y,n,t),  the  body  force  per  unit  volume  by 
b(y,t).  and  the  mass  per  unit  volume  by  p{y,i).  Also,  the  Eulerian  form  of  the 
velocity  of  the  particle  of  material  at  y  e  Rt  at  time  f  e  P  is  denoted  by  v(y ,  t).  The 
familiar  Eulerian  global  form  of  the  balance  of  mass,  linear  momentum,  and  angular 
momentum  are 


(2.1) 


J  y  X  tdA  +  J  y  X  hdV  =  -^  J  y  x  p^dV, 


i22) 


(23) 


dDf  Dt  Dt 

respectively,  V  Dt  C  Rt  and  V  t  €  P.^  Equations  (2.1),  (2.2),  and  (2.3)  can  be  used 
to  derive  the  rate  of  woik-energy  equation  given  by 


*  In  this  paper,  math-italk  quantities  denote  scalan  and  bold-faced  quantities  denote  lenson,  including  vectors. 
^  In  this  paper,  unless  otherwise  stated,  whenever  a  subset  of  a  tegular  region  is  considered,  it  is  assumed  that 
the  subset  is  a  regular  subregioo.  Similaiiy,  all  subregions  of  regular  regions  are  assumed  to  be  regular. 
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J  t  •  vdA  +  y  b  •  \dV  =  jr^  VyXdV  +  ^  ^  ipv  •  \dV,  (2.4) 

l^t  Dt  Df 

V  Dt  C  Rt  and  V  t  e  F,  where  T(y,<)  is  the  true  (or  Cauchy)  stress  tensor  and  VyV 
denotes  the  gradient  of  9(y,  t)  with  respect  to  y. 

3.  Multiple  Reference  Configurations  and  the  Continuum  Model 

Consider  a  region  R  in  E3  that  the  body  B  can  occupy,  in  the  sense  that  there 
exists  a  suitably  smooth  and  inveitible  mapping  that  maps  R  into  Rt.  Note  that  R  is 
not  such  that  the  body  B  has  to  occupy  it  at  some  time  t  £  F.  Such  a  region  R  can 
be  used  as  a  reference  configuration  for  the  body  B. 

In  the  continuum  mode!  that  is  developed  in  this  paper,  each  phase  has  its 
own  constitutive  relation.  These  constitutive  relations,  however,  are  not  all  defined 
with  respect  to  the  same  fixed  reference  configuration.  Instead,  each  phase  has  its 
own  reference  configuration  for  the  definition  of  its  constitutive  equations,  and  for 
the  expression  of  its  field  equations.  More  specifically,  for  the  continuum  model 
that  is  develped,  it  is  assumed  that  each  phase  has  a  configuration  corresponding 
to  an  unstressed  undeformed  state,  and  that  each  phase  behaves  elastically  for 
some  range  of  deformations  about  its  unstressed  undeformed  configuration.  The 
reference  configuration  for  each  phase  of  these  materials  is  taken  to  coincide  with  the 
unstressed  undeformed  configuration  of  that  phase.  With  respect  to  the  undeformed 
austenite  phase,  these  reference  configurations  coincide  the  shape  deformations  of 
the  martensites. 

For  simplicity,  in  the  following  sections,  the  continuum  model  will  be  developed 
and  the  corresponding  field  equations  and  jump  conditions  will  be  derived  for  a 


twO'phase  material.  The  corresponding  results  for  a  material  that  consists  of  mote 
than  two  phases  can  be  obtained  in  a  similar  manner. 

4.  The  Kinematics  Using  Multiple  Reference  Configurations 

Consider  the  body  B  that  was  described  in  Section  2.  Additionally,  assume  that 
this  body  consists  of  two  phases,  which  will  henceforth  be  referred  to  as  phase  1  and 
phase  2.  Assume  that  at  each  t  €  P  phase  1  occupies  a  subregion  of  Rt,  and 
phase  2  occupies  a  subregion  R^  of  R*,  where  Rf  U  R^  =  Rt  and  Rf  n  R^  =  0 
(Hgure  1).  These  two  subregions  of  Rt  are  separated  by  an  interface  St  which  can 
pass  over  particles  of  material  in  Rt.  If  this  occurs,  Rf  will  increase  in  size  while 
R^  decreases  in  size,  or  vise-versa,  depending  on  the  direction  of  motion  of  St. 

It  is  assumed  that  there  exists  a  configuration  R  of  phase  1  that  corresponds  to 
an  unstressed  undeformed  state  of  that  phase.  This  configuration  R  will  be  used  as 
a  stationary  reference  configuration  for  Rt  (Hgure  1).  Let  x  denote  a  point  or  the 
position  vector  of  a  point  in  R.  Let  y(x,t)  be  the  suitably  smooth  and  invertible 
mapping  which  maps  R  into  Rt  at  each  t  G  P,  with  y  =  y(x,  t)  =  x+u(x,  t)  V  (x,t)  € 
Rx  P,  where  u  =  u(x,  t)  is  the  displacement  of  the  point  y  =  y(x,  t)  firom  the  point 
X  at  time  t  €  P.  The  deformation  gradient  of  y  is  defined  as  F(x,f)  =  Vy(x,t), 
and  the  Jacobian  of  F  is  defined  as  J(x,t)  =  detF(x,t),  where  it  is  required  that  J 
>  0  to  exclude  reflections.^  The  velocity  of  the  particle  of  material  at  y  =  y(x,  t)  is 
defined  as  v(x,  t)  =  ^y(x,  t),  and  the  Eulerian  form  of  the  velocity  used  in  Section 
2  can  be  defined  as  v(y,t)  =  v(x(y,f),t),  where  x(*,t):Rt  R  at  each  f  €  P 
is  the  inverse  of  y.  Let  R“  =  x(Rp,<),  R'*'  =  x(R^,<),  and  S  =  x(St,<).  Note 
that  R~  U  R'*'  =  R,  R~  n  R"*"  =s  0,  S  is  the  surface  separating  R~  firom  R'*’,  and  S 
moves  within  R  as  St  moves  within  Rt. 

*  V  denotes  the  gradient  operator  with  respect  to  x  €  R. 


It  is  assumed  that  there  exists  a  shape  deformation  for  phase  2  with  respect 
to  R  which  corresponds  to  an  unstressed  undeformed  configuration  of  that  phase. 
Let  Rf  be  the  reference  configuration  coinciding  with  this  sh^  deformation  of 
phase  2  for  all  t  €  r  (Figure  1).  We  note  that  because  R  has  been  assumed  to  be 
stationary  and  there  must  be  continuity  of  displacements  at  the  phase  boundary  for 
all  t  e  F,  R^  will  most  likely  be  moving  if  a  problem  other  than  a  static  problem 
is  considered.^  Let  Xi  denote  a  point  or  the  position  vector  of  a  point  in  R^,  and 
let  Xi(x,t)  be  the  suitably  smooth  and  invertible  mapping  that  maps  R*^  into  Rf  at 
each  t  €  F,  with  Xi  =  Xi(x,t)  VxgR'*'  at  each  i  €  F.  This  moping  xi  is  assumed 
to  be  given  in  a  problem.  Let  F(x,t)  =  Vxi(x,t)  and  /(x,<)  =  detF(x,<),  with 
J  >  0.  Let  yi(xi,t)  be  the  suitably  smooth  and  invertible  mapping  that  maps 
Ri  into  R^  at  each  <  6  F,  with  y  =  yi(xi,i)  =  X|  +  Ui(xi,t)  V  Xi  €  Rj" 
at  each  <  €  F.  Note  also  that  y  =  yi(xi(x,t),t)  =  y(x,<)  Vx  G  R"*"  at  each 
t  G  F.  Let  Fi(xi,t)  =  Viyi(xi,<)  and  Ji(xi,t)  =  detFi(xi,t),  with  Ji  >  0.* 
The  velocity  field  of  the  particles  of  material  in  phase  2  as  a  function  of  xi  G  Rf 
is  given  by  Vi(xi,  t)  =  v(yi(xi,  <),  t)  V  XiG  Rf  at  each  t  G  F,  or  equivalently  as 
=  v(x(xi,t),<)  yxiGRi"  at  each  f  G  F,  where  x(-,t):Ri'  •-♦R'''  at  each 
t  G  F  is  the  inverse  of  Xi.  Also,  Xi(',t):R^  R^  at  each  f  G  F  represents  the 
inverse  of  yi. 

Let  N(x,t)  represent  a  unit  vector  normal  to  S  that  points  into  R'*',  and  let 
L(x,  t)  represent  a  vector  tangent  to  S,  both  at  a  point  on  S  coinciding  with  the  point 
X  G  R  at  time  <  G  F.  Also,  let  V(x,<)  represent  the  velocity  of  the  point  on  S 
coinciding  with  the  point  x  G  R  at  time  t  G  F.^  In  the  following,  if  ^(x,  t)  represents 

^  See  Section  6.  Equation  (6.9).  for  an  example  of  this. 

'  Vi  denotes  the  gradient  opentor  with  respect  to  xt  6  R|’* 

*  In  the  rest  of  this  thesis,  a  point  on  the  surface  S  that  coincides  with  the  point  x  C  R  at  lime  (  £  T  will  simply 
be  referred  to  as  the  point  x  €  S.  Also,  note  that  V  represents  a  nominal-type  velocity  and  not  the  velocity  of  a 
point  on  St. 
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a  generic  field  quantity  that  is  discontinuous  at  S,  then  and  g'^(x,t)  denote 

the  limiting  values  of  ^  at  x  €  S  as  this  point  x  is  approached  from  negative  and 
positive  sides  of  S,  respectively,  in  directions  parallel  to  N(x,t).‘° 

Because  the  displacements  at  the  interface  separating  two  phases  involved  in  a 
martensitic  phase  transformation  are  continuous  while  the  strains  at  the  interface  may 
be  discontinuous,  we  require  that  y(x,t)  be  continuous  on  R  x  P,  and  allow  the  first 
and  second  derivatives  of  y  to  be  piecewise  continuous  on  R  x  P,  with  discontinuities 
occurring  only  at  points  on  S.  As  a  result  of  this,  we  have 

y+(xi(x,  t),  t)  =  y~(x,  <),  (4.1) 

Vx  6  S  at  each  f  6  P.  Taking  the  differential  of  both  sides  of  (4.1)  while  keeping 
time  fixed  (and  recalling  the  continuity  conditions  on  y)  yields 

(f+F+-F-)l  =  0,  (4.2) 

V  X  €  S  at  each  t  G  P,  and  for  every  vector  L  tangent  to  S  at  x  e  S,  where 
Xi  =  xi(x,t)  in  F^.  Differentiating  (4.1)  with  respect  to  time  yields 

v+  -  V-  +  (f+F+  -  F-)v  =  0,  (4.3) 

Vx  €  S  at  each  f  G  P,  where  Xi  =  Xi(x,f)  in  and  FJ".” 

It  can  be  shown  that  if  given  an  F'  ,F^,  and  an  F'*'  such  that  Equation  (4.2) 
is  satisfied  at  a  point  on  some  surface,  there  exists  vectors  a  and  N  defined  at  that 
point  such  that 

F+F+-F"=a®N,  (4.4) 

If  ■  pieoewue  condnuous  6eld  quantity  is  discontinuous  at  S.  it  is  sometimes  said  that  the  quantity  jumps 
across  S,  and  an  equation  relating  g~  and  g*  is  sometimes  referred  to  u  a  jump  condition. 

"  The  conesponding  forms  of  Equations  (4.2)  and  (4.3)  involving  only  one  reference  configuration  are  well 
known  and  can  be  found  in.  e.g..  [2]. 
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where  N  is  normal  to  the  surface  at  that  point.'^  Also,  note  that  a  0  N  in  this 
equation  is  a  rank-one  two-tensor. 

5.  The  Nominal  Form  of  the  Field  Equations  Using  Multiple  Reference 
Configurations 

In  addition  to  the  continuity  requirements  on  y  and  u  discussed  in  the  previous 
section,  it  is  assumed  that  b  is  continuous  V  y  €  Ri  and  V  t  G  F,  and  that  r  and 
its  gradient  are  piecewise  continuous  V  y  €  Rt  and  V  t  €  F,  with  discontinuities 
occurring  only  at  points  on  Sf 

The  global  form  of  the  balance  of  mass  given  by  Equation  (2.1)  can  be  expressed 
with  respect  to  R  as 

i/^  =  0.  (5‘) 

D 

where  D  =  x(Dt,t),  and  p  —  Jp  represents  the  mass  per  unit  volume  of  R.  Note 
that  since  (2.1)  is  valid  VDt  C  Rt,  (5.1)  is  valid  V  D  c  R.  We  require  that  p  be 
continuous  Vx  G  R.  Localizing  Equation  (5.1)  yields  the  familiar  result  that  p  must 
be  independent  of  time  V  x  G  R.  Thus, 

p(x)  =  J(x,t)^y(x,t),0,  (5^) 

V(x,  i)  G  R  X  F.  Note  that  the  global  field  equations  given  by  (2.1)-(2.3)  are  still  valid 

in  the  regions  of  space  and  time  indicated  there,  for  the  case  where  p,  r,  and  b  have 

the  continuity  conditions  specified  in  this  section  and  y  has  the  continuity  conditions 

specified  in  the  previous  section.  However,  for  these  continuity  conditions,  the  wmk- 

energy  equation  given  by  (2.4)  is  valid  only  for  subregions  of  Rt  not  containing  a 

portion  of  St,  as  will  be  discussed  further  in  Section  8. 

Hie  conespondins  form  of  this  equation  with  F^F'*’  replaced  with  T*,  where  F'*’  =  FfF'*’,  is  well  known 
and  can  be  found  in,  e.g..  [14]. 
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For  subregions  of  R^,  the  global  forms  of  the  balance  of  linear  and  angular 
momentum  for  phase  1  given  by  (2.2)  and  (2.3),  respectively,  can  be  expressed  with 
respect  to  F.  as 


I 


3D- 


andA  + 


/ 


D- 


J 


ao- 


y 


andA  A-  J  y  fdV  = 

D- 


(5.3) 


(5.4) 


respectively,  where  D“  =  x(Df,<),  <r(x,t)  =  /(x,t)r(y(x,<),0F"^(x,<)  is  the 
nominal  stress  tensor  with  respect  to  R“,  f(x,t)  =  J(x,<)b(y(x,t),t)  is  the  nominal 
body  force  per  unit  volume  of  R~,  and  n(x)  is  the  outward  unit  normal  vector  field 
on  the  boundary  of  D“.  These  nominal  field  equations  can  be  obtained  in  the  usual 
way  from  Equations  (2.2)  and  (2.3),  respectively  (see,  e.g.,  [13]). 

For  subregions  of  R^^,  the  global  forms  of  the  balance  of  linear  and  angular 
momentum  for  phase  2  given  by  (2.2)  and  (2.3),  respectively,  can  be  expressed  with 
respect  to  R]^  as 

(5.5) 

(5.6) 


J  axnxdA+  j  fxdV  = 


8D+ 


D* 


J  Si  X  <rin,dA+  j  St  ^  ^t^  -  ^  Pi^t^V^ 

an*  D* 


respectively,  where  D|  =  Xi(D^,t),  Pi(xj,<)  =  Ji(xift)p(yi(xi,i),t)  is  the  mass 
per  unit  volume  of  RJ",  o'i(xi,t)  =  Ji(xi,t)Ti(yi(xi,t),<)F7^(X|,t)  represents  the 
stress  tensor  with  respect  to  R^,  fi(X|,f)  =  Ji(xi,t')b(yi(xi,i),t)  represents  the 
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body  force  per  unit  volume  of  R|,  and  ni{xi,i)  is  the  outward  unit  normal  vector 
field  on  the  boundary  of  D^.  Also,  pi  is  related  to  p  by 

p(x)  =  J(x,t)^i(xi(x,t),t),  (5.7) 


V  X  €  R'*’  at  each  t  €  F.  Note  that  a-y,  f„  and  pi  do  not  have  a  true  nominal  form 
since  Xi  is  a  function  of  time.  In  fact,  because  of  this  time  dependence  these  quantities 
are  closer  to  having  an  Eulerian  form.  Note  also  that  Equations  (5.5)  and  (5.6)  can 
be  obtained  from  Equations  (2.2)  and  (2.3)  in  a  manner  completely  analogous  to  that 
used  to  obtain  Equations  (5.3)  and  (5.4),  regardless  of  whether  Rf  is  stationary  or 
whether  points  in  R^  are  moving,  which  is  the  case  if  x^  is  a  function  of  time. 

Localization  of  the  global  nominal  balance  laws  given  by  (5.3)  and  (5.4)  using 
(5.2)  yields 


diver  4-  f  =  pa, 

=  Fa'^, 


(5.8)» 


respectively,  Vx  €  R~  at  each  <  €  F,  where  a(x,<)  =  ^v(x,<)  is  the  acceleration 
of  the  particle  of  material  in  RT  that  corresponds  to  the  point  x  e  R~. 

Localization  of  the  global  nominal  balance  laws  given  by  (5.5)  and  (5.6)  using 
(5.7)  yields 


divi<7-i  +  fi  =  piai, 


(5.9)*^ 


respectively,  Vxi  €  Rj"  at  each  t  €  F,  where  ai(xi,t)  =  [^Vi(xi(x,t),t)]j(j^^ 


div  denotes  the  divergence  opemtor  with  respect  to  x  €  R. 
divi  denotes  the  diveigenoe  operator  with  respect  to  xi  c  R^. 


12 


is  the  acceleration  of  the  particle  of  material  in  that  corresponds  to  the  point 
Xi  =  Xi(x,t)  in  R^. 

For  subregions  Dt  of  Rt  which  contain  a  portion  of  St,  the  global  form  of  the 
balance  of  linear  momentum  given  by  Equation  (2.2)  can  be  expressed  with  respect 
to  R  in  terms  of  the  field  quantities  defined  on  R~  and  R^.  Localizing  this  equation 
at  points  on  S  would  then  yield 

-  (T-^N  +  p(v+  -  v-)(V  .  N)  =  0,  (5.10) 


Vx  e  S  at  each  t  €  T,  where  Xi  =  Xi(x,<)  in  and 

The  balance  of  angular  momentum  given  by  Equation  (2.3)  for  a  subregion  Dt 
of  Rt  containing  a  portion  of  St  is  automatically  satisfied  if  the  jump  conditions 

(4.2) ,  (4.3),  and  (5.10)  are  satisfied  at  the  points  on  S  corresponding  to  the  points 
on  the  portion  of  St  contained  in  Dt.  Also,  the  local  field  equations  given  by 

(5.2)  and  (5.7)-(5.9),  the  kinematic  jump  condition  (4.1),  or  (4.2)  and  (4.3),  and  the 
linear  momentum  jump  condition  (5.10)  are  all  together  equivalent  to  the  global  field 
equations  given  by  (2.1)-(2.3). 

6.  Elastic  Materials  and  Multiple  Reference  Configurations 

In  the  following,  it  is  assumed  that  phase  1  behaves  elasticaUy  for  some  range 
of  deformations  about  R~  and  that  phase  2  behaves  elastically  for  some  range  of 
deformations  about  R|.  In  particular,  it  is  assumed  that  phase  1  possesses  an  elastic 

”  Tile  corretpondinf  form  of  this  jump  condition  invoiving  only  one  reference  configuration  it  well  known 
and  can  be  found  in.  e.g.,  (2].  In  fact,  the  jump  conditions  (4.2),  (4.3),  and  (S.IO)  can  formally  be  obtained 
from  the  jump  conditiofu  presented  in  [2]  by  replacing  F'*',  v'^,  and  with  9*,  and  , 

respectively. 


13 


potential 

W^  =  W^(P(x,<),x),  (6.1) 

Vx  €  R~  and  VF  €  such  that  the  nominal  stress  tensor  for  this  phase  is  given  by 

a(x,  t)  =  Wr(F(x,  <),  x),  (6.2) 

Vx  e  R“  and  VF  e  where  £■*■  is  a  subset  of  C^,  is  the  set  of  all  two-tensors 
with  positive  determinants,  and  represents  the  range  of  deformations  about  R~ 
for  which  phase  1  behaves  elastically.  Equations  (6.1)  and  (6.2)  correspond  to  the 
standard  definition  of  an  elastic  (or  hyperelastic)  material  defined  with  respect  to  a 
single  fixed  reference  configuration,  that  can  be  found  in  almost  any  textbook  on 
finite  elasticity  (see,  e.g.,  [13]).  The  constitutive  equations  (6.1)  and  (6.2)  can  have 
these  standard  forms  because  the  reference  configuration  R"  is  stationary.  In  fact, 
as  a  result  of  R~  being  stationary  and  (6.1)  and  (6.2),  the  general  field  equations  in 
terms  of  the  displacements  for  phase  1,  and  consequently  the  linearized  constitutive 
equations  and  the  linearized  field  equations  in  terms  of  the  displacements  for  phase 
1,  will  all  have  the  standard  well  known  forms  corresponding  to  those  for  elastic 
materials  where  single  fixed  reference  configurations  are  used. 

For  phase  2,  we  assume  that  there  exists  an  elastic  potential  Wi  defined  with 
respect  to  RJ*'.  It  is  assumed  that  is  a  function  of  F,.  Additionally,  as  one 
might  expect,  it  is  required  that  the  inhomogeneity  of  Wi  remain  the  same  for  each 
particle  of  material  in  phase  2  as  time  progresses.  If  xi  is  a  function  of  time,  the 
reference  point  Xi  €  R^  for  a  given  particle  of  material  in  phase  2  is  changing  as 
time  progresses,  and  a  different  particle  of  material  occupies  a  given  point  Xi  €  R| 
at  different  times  t  €  F.  Because  of  this,  the  inhomogeneity  of  Wi  cannot  be 
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expressed  with  respect  to  points  xi  €  Rt  if  xi  is  a  function  of  time.^^  Instead,  the 
inhomogeneity  of  Wi  must  be  represented  with  respect  to  points  x  in  the  stationary 
reference  configuration  R'*',  so  that  the  effect  of  the  inhomogeneity  of  Wi  on  a  given 
particle  of  material  in  phase  2  follows  that  particle  for  all  t  G  F.  Thus,  we  assume 
that  Wi  has  the  form 


W^i  =  W^t(Fi(x„t),x(x„t)),  (6.3) 

V  Xi  e  Ri",  V  #  6  r,  and  V  Fi  €  ,  where  is  a  subset  of  and  represents 

the  range  of  deformations  about  R|  for  which  phase  2  behaves  elastically.  It  is 
further  assumed  that  Wi  is  such  that  the  stress  tensor  for  phase  2  with  respect  to 
R^  is  given  by 


<Ti(xi,t)  =  H^ij.^(Fi(xi,t),x(xi,t))  (6.4) 

Vxi  €  R^,  Vt  €  r,  and  VFi  €  JCf. 

As  discussed  in  the  Introduction,  the  symmetry  of  the  crystal  lattice  of  a  material 
is  represented  in  the  material  symmetry  group  of  that  material  (see  [13]).  It  is  required 
that  the  elements  of  a  material  symmetry  group  of  a  solid  be  unimodular,  so  that  they 
preserve  volume.*’  The  material  symmetry  group  for  phase  1  is  defined  with  respect 
to  R~,  and  the  material  symmetry  group  for  phase  2  is  defined  with  respect  to  Rf. 
These  material  symmetry  groups  restrict  tl»  functional  forms  of  the  elastic  potentials 
of  their  respective  phases.  Mote  specifically,  if  ^(x)  is  the  material  symmetry  group 
for  phase  1,  IF  is  required  to  be  such  that 


IF(FH,x)  =  IF(F,x),  (6J) 


However,  if  xi  is  not  i  function  of  tune,  which  would  be  appropriate  for  the  static  case,  the  referenoe 
configuration  Rf  is  stationary,  and  consequently  the  inhomogeneity  of  Wi  can  be  expressed  explicitly  with 
respect  to  points  xi  €  Rf- 

By  definition,  a  unimodular  two-tensor  has  a  determinant  equal  to  one. 
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V  H  €  0(x)  and  V  F  €  at  each  x  €  R“.  Similarly,  if  ^i(x)  is  the  material 
symmetry  group  for  phase  2,  W\  is  required  to  be  such  that 

M^i(F,H|,x)  =  W^,(Fi,x),  (6.6) 

V  Hi  €  ^i(x)  and  V  Fi  €  £■*■  at  each  x  e  R+.  We  note  that  because  a  different 
elastic  potential  is  used  for  each  phase,  the  material  symmetry  groups  for  phase  1 
and  phase  2  can  be  chosen  independently  of  each  other  and  aihitrarily,  with  the 
exception  that  they  must  be  subsets  of  the  unimodular  group.  Consequently,  the 
material  syrrunetry  group  of  each  phase  can  be  chosen  to  reflect  any  type  of  crystal 
symmetry  with  any  orientation.  Therefore,  the  change  in  crystal  structure  that  takes 
place  during  a  martensitic  phase  transformation  and  the  orientation  relation  between 
the  crystal  lattices  of  the  austenite  and  the  martensite  can  be  accurately  represented. 
In  fact,  because  the  elastic  potential  for  each  phase  is  defined  with  respect  to  the 
undeformed  configuration  of  that  phase,  the  material  symmetry  group  for  that  elastic 
potential  can  be  chosen  to  be  the  crystallographic  point  group  corresponding  to 
the  crystal  symmetry  of  that  phase.  This  certainly  makes  constructing  the  elastic 
potential  for  each  phase  much  more  feasible. 

We  next  require  that  Xi  has  the  form 

Xi(x,f)  =  ^(x)  +  0(<),  (6.7) 

V  X  €  R***  at  each  f  €  F.  If  Xi  has  this  form,  the  constitutive  behavior  for  phase 
2  given  by  (6.3)  and  (6.4)  results  in  the  stress-power  for  any  subregion  of  phase  2 
being  equal  to  the  time  rate  of  change  of  the  integral  of  Wi  over  that  subregion;  i.e. 
if  Xt  has  the  form  given  by  (6.7), 

J  <ri .  =  WrdV, 

D*  D* 


(6.8) 
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VD^  c  Rj*  and  Vr  e  F  (see  [9]).  For  this  case,  it  can  easily  be  shown  that  the  change 
in  the  total  strain  energy  and  kinetic  energy  of  any  subregion  of  phase  2  in  a  complete 
cycle  is  zero.  This  is  certainly  the  most  important  property  of  an  elastic  material,  and 
any  choice  of  Xj  and  the  forms  of  and  the  corresponding  stress  tensor  a,  that  did 

not  result  in  this  would  not  be  consistent  with  the  standard  definition  of  an  elastic 
material.'*  For  the  important  special  case  where  x,  is  a  homogeneous  deformation 

and  there  is  continuity  of  displacements  at  the  phase  boundary  when  F  =  lVxeR’ 
and  Fj  =  IVx^  €  R*,  at  each  t  e  F,  x,  must  have  the  form 

x,(x,r)  =  x+([x-x(r)]n)a,  (6.9) 

where  x(r)  is  any  point  on  the  interface  S,  and  a  and  n  are  constant  vectors.  We 
note  that  Xj  in  (6.9)  has  the  form  of  an  invariant  plane  strain  [15].  Additionally,  for 
the  Xj  given  by  (6.9)  and  when  F  =  IVx  €  R'  and  Fj  =  IVXj  e  R^,  at  each  r  e  F,  the 
interface  S  is  a  plane  with  unit  normal  n  and  translates  with  velocity  V(t)  =  ^x(r), 
and  all  points  in  phase  2  translate  with  velocity  v(r)  =  -(V(f)  0)3. 

7.  The  Domains  of  the  Elastic  Potentials 

As  indicated  in  the  previous  section,  the  elastic  potentials  given  by  (6.1)  and  (6.3) 
are  defined  for  finite  deformations  in  2^  czL*  and  respectively. 

Additionally,  2^  is  with  respect  to  R,*  in  the  sense  that  the  two-tensors  in  2^ 

Note  also  that  if  X^  has  the  form  given  by  (6.7),  the  transformation  strain  for  phase  2  is  independent 
of  time,  which  is  probably  most  appropriate  for  the  case  considered  here  (i.e.,  for  the  case  where 
phase  2  represents  a  phase  that  is  in  a  purely  mechanical  process  and  behaves  elastically). 
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represent  deformations  Fi  with  respect  to  The  corresponding  subset  of  that 
represents  with  respect  to  R  can  be  defined  as 

£t  =  {f/  F  =  FiP,  where F,  e  }.  (7.1) 


Let  S*  denote  the  set  of  all  symmetric  positive  definite  two-tensors.  In  the 
space  5'*',  the  deformations  with  respect  to  R  corresponding  to  the  undeformed 
configurations  of  phase  1  and  phase  2  are  C  =  1  and  C  =  F^F,  respectively. 
Considering  this,  the  sets 


5+  =  |c/  C  =  F'^F,  where  F  €  £"■}, 
St  =  {C/C  =  FF,  where F  6  £}■ } 


(7.2) 


represent  the  domains  with  respect  to  R  about  0  =  1  and  C  =:  F'^F  in  S'*'  that 
phase  1  and  phase  2,  respectively,  behave  elastically  for.  We  note  that  if  S'*'  and 
St  are  relatively  large  and  1  and  C  are  relatively  close,  these  two  phase  domains 
may  intersect  in  S'*'.  In  this  case,  some  might  consider  this  constitutive  model  to  be 
multivalued.  Either  way,  an  intersection  of  these  two  phase  domains  will  not  present 
any  difficulties  in  setting  up  or  solving  a  boundary  value  problem. 


Most  metals,  however,  have  phases  that  are  not  elastic  for  finite  deformations 
about  their  undeformed  configurations.  They  instead  have  phases  that  have  elastic- 
plastic  constitutive  behaviors  with  yield  stresses  corresponding  to  infinitesimal 
deformations.  If  the  material  under  consideration  represents  such  a  material  and 
if  the  stresses  within  phase  1  and  phase  2  are  less  than  the  yield  stresses  of  these 
phases,  respectively,  for  all  f  €  F,  both  of  these  phases  will  behave  elastically  in  this 
time  interval.  Additionally,  in  this  case,  because  the  deformations  are  infinitesimal 
for  all  f  €  F,  the  linearized  forms  of  the  constitutive  equations  given  by  (6.1) 
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and  (6.2)  about  F  =  1  and  (6.3)  and  (6.4)  about  Fi  =  1  can  be  used.  This 
is  probably  most  appropriate  in  the  temperature  interval  containing  the  martensitic 
start  temperature,  where  the  stress  necessary  to  induce  nucleation  or  growth  of  a 
variant  of  martensite  is  less  than  the  yield  stresses  of  both  phases.  In  this  case, 
the  material  will  deform  by  a  martensitic  phase  transformation  before  it  wiU  deform 
by  plastic  deformation.^^  This  martensitic  phase  transformation  may  be  associated 
with  the  austenite  phase  transfonning  into  a  martensite  phase,  or  with  one  variant 
of  martensite  transforming  into  another  variant  of  the  same  martensite,  which  is 
also  known  as  reorientation.  Also,  for  this  case  where  and  Cf  contain  only 
deformations  such  that  IVu]  «  1  and  |ViUi|  «  1,  respectively,  if  Xi  is  a  finite 
deformation,  there  is  no  chance  that  S'*'  and  5j  will  intersect 


8.  The  Driving  Traction 

As  mentioned  in  Section  5,  the  global  form  of  the  rate  of  work-energy  equation 
given  by  (2.4)  is  not  valid  for  subregions  of  Rt  containing  a  portion  of  St  since  the 
continuity  assumptions  that  are  necessary  for  the  derivation  of  that  equation  from 
Equations  (2.1)-(2.3)  do  not  exist  for  these  subregions.  For  the  two-phase  elastic 
material  under  consideration,  the  global  form  of  the  rate  of  work-energy  equation  for 
subregions  Dt  of  Rt  containing  a  portion  St  of  St  can  be  expressed  as 


/' 


t  •  vdA  -h  /  b  •  xdV 

Dt 


=  - / 
dt  J 


D- 


WdV  +  -^  J  WidV  +  ^Jpy-  vdv, 

Dt 


(8.1) 


A  maieriil  defoniu  by  a  manensitic  phase  transformation  in  the  sense  that  the  creatioo  of  the  manensite  phase 
produces  a  deformation  due  to  the  shape  deformation  of  that  martensite. 
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where  D"  =  x(D^,i),  Df  =  S*  =  x(St*,<),  U  D?"  =  Dj,  and 


/  =  (j+Ty+  -  W-')  -  +  a- j  •  (f+F+  -  F-). 


(8^) 


The  quantity  /N  is  referred  to  as  the  driving  traction,  and  /  as  the  scalar  driving 
traction.^  The  integral 


(8.3) 


can  be  interpreted  as  representing  the  rate  of  work  done  on  the  interface  S*  by  the 
traction  /N  exerted  by  the  body  on  the  interface.  We  require  that  this  integral  be 
positive  so  that  it  represents  a  dissipation  of  energy.  Localizing  this  equation  at 
points  on  the  interface  S  then  yields 

/N.V>0  (8.4) 

V  X  e  S  at  each  t  e  Note  that  energy  can  be  dissipated  only  at  the  phase 
boundary,  and  that  if  /N  -  V  =  0VxeSat  each  <  €  F,  energy  is  conserved.  Thus, 
by  allowing  /N  •  V  >  0  we  are  in  effect  considering  a  nonconservative  system  even 
though  both  phases  behave  elastically  for  points  not  on  the  interface. 

We  can  also  postulate  a  constitutive  relation  relating  the  scalar  driving  traction 
to  the  normal  component  of  the  phase  boundary  velocity  at  each  x  G  S  (see  [2]). 

^  Refer  to  [2]  and  [S]  for  more  extensive  discussions  about  the  driving  tractioo.  Also,  the  scalar  driving 
tractioo  given  by  (8.2)  can  formally  be  obtained  from  the  scalar  driving  tractioo  derived  in  [2]  by  replacing 
W*,  a*,  and  with  and  PJ’F'*’,  respectively. 

This  requirement  is  also  equivalent  to  the  second  law  of  thermodynamics  for  the  type  oi  prooeu  under 
consideratioo  (see  [2]). 
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More  specifically,  we  can  postulate  that 

V.  =  $(/,x),  (8^) 

at  each  x  €  S,  where  $(/,  x)  is  given  and  depends  on  the  material,  and  Vn  =  V  •  N. 
This  may  be  done  primarily  to  provide  an  extra  equation  for  the  extra  unknowns  — 
the  variables  describing  the  location  of  the  phase  boundary. 

9.  The  Linearized  Problem 

As  mentioned  in  the  Introduction,  one  of  the  main  advantages  of  using  the 
continuum  model  that  was  developed  in  this  paper  is  that  the  field  quantities  and 
equations  are  in  forms  that  permit  direct  linearization,  while  retaining  finite  shape 
deformations  for  the  martensite  phases.  This  is  the  case  since  the  displacements 
for  each  phase  are  measured  from  the  reference  configuration  coinciding  with  an 
unstressed  undeformed  configuration  of  that  phase  (i.e.  the  shape  deformation  of  that 
phase),  and  consequently,  for  the  appropriate  boundary  and  initial  conditions,  the 
displacement  gradients  can  be  considered  infinitesimal.^  Another  advantage  of  using 
this  continuum  model  is  that  for  the  linear  case  the  nominal  stress  for  each  phase 
is  approximately  equal  to  the  true  stress  for  that  phase.  This  is  very  convenient  for 
solving  certain  types  of  boundary  value  problems,  as  will  become  more  apparent  in 
Sections  13  and  14  where  several  of  these  types  of  problems  are  considered. 

In  the  next  few  sections,  the  linearized  field  equations  and  jump  conditions  in 
terms  of  the  displacements  for  the  two-phase  elastic  material  under  consideration  are 
considered.  The  general  equations  in  terms  of  the  unspecified  elastic  potentials  given 
by  (6.1)  and  (6.2)  for  finite  deformations  can  be  obtained  in  a  similar  manner  (see  [9]). 

”  b  is  also  assumed  here  that  the  unstressed  undefonned  configuiaiioo  of  each  phase  cofrespoods  to  a  local 
minimum  of  the  elastic  potential  for  that  phase. 
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10.  The  Linearized  Constitutive  Equations  and  Field  Equations 

In  the  following,  it  is  assumed  that  the  undefonned  configurations  of  phase  1 
and  phase  2  correspond  to  local  minima  of  their  respective  elastic  potentials.  For  the 
linear  problem,  we  assume  that  |Vu|  «  1  Vx  6  R“  and  |Vui|  «  1  V Xi  €  Rf  at 
each  t  e  r.  In  this  case,  we  have  for  phase  1 

IV  =  ly-  +  iVu  •  (CVu)  +  0(|Vuf) 

(10.1) 

=  lV  +  ic-(Ce)  +  0(|c|®) 

ff  =  CVu  +  0(|Vup) 

=  Ce  +  0(|cp) 

^(y(x,t),t)  =  <r(x,f)  +  O(lVup) 

<T  =  <rT^.o(|Vu|=^) 

Vx  €  R~  at  each  t  eT,  where  e  =  5  ^Vu  +  (Vu)^^  is  the  infinitesimal  strain  tensor 
for  phase  1.  lV*(x)  =  lV(l,x),  and  C(x)  with  C,j«  =  |  is  the  elasticity 

four-tensor  for  phase  1.  For  phase  2,  we  have 

=  w:  +  IV, u^ .  (C,V,Ut)  +  0(|V,u,n 

(10.2) 

=  lV:-HiMC,ei)  +  0(|ei  ft 

<ri  =  C,V,u,  +  0(|V,u,|*) 

=  Cifij  -f-  0(1  Cl  p) 

’•i(yi(x,,<),<)  s=  o’i(xi,t)  4-  0(|ViUi|*) 


22 


^l=4  +  0(|VtUxft 


Vxi  €  Ri"  at  each  t  G  F,  where  e\  =  is  the  infinitesinial  strain 

tensor  for  phase  2,  W^f(x(xi,t))  =  1^1(1, x(xi,t)),  and  Ci(x(xi,t))  with  = 
is  the  elasticity  four-tensor  for  phase  2.  The  elasticity  tensors  C 


Fi=l 


and  Cl  are  both  four-tensors  which  contain  the  material  coefficients  and  reflect  the 
material  symmetry  of  phase  1  and  phase  2,  respectively.  Additionally,  both  C  and 
Cl  are  positive-definite,  since  it  has  been  assumed  that  the  undeformed  configuration 
of  each  phase  corresponds  to  a  relative  minimum  of  the  elastic  potential  for  that 
phase.  These  four-tensors  also  have  the  usual  major  and  minor  symmetries. 


Substituting  the  constitutive  equation  (10.1)2  into  the  equation  representing  the 
balance  of  linear  momentum  for  phase  1  given  by  (5.8)i  and  expressing  a  in  that 
equation  in  terms  of  u  yields,  in  a  given  coordinate  frame. 


dCjjki  duk  ,  ^  ^ 

^  Q  I  Q  «  I  Ji  P  ’ 
oxj  ox\  oxioxj  at^ 


(10.3) 


Vx  G  R"  at  each  t  G  F,  where  the  second-order  terms  have  been  neglected. 

Similarly,  substituting  the  constitutive  equation  (10.2)2  into  the  equation 
representing  the  balance  of  linear  momentum  for  phase  2  given  by  (5.9)i  and 
expressing  ai  in  that  equation  in  terms  of  xi  and  Ui  yields 


divi(CiViUi)  -I-  fi  =  (xi  -f-  ui)  , 


(10.4) 


Vxi  G  Rf  at  each  t  G  F,  where  Xi  =  ^Xi(x,t)  =  ^©(f),  Uj  = 
I  ctUi(xi(x,  <),<)(  ,  and  the  second-order  terms  have  been  neglected.  In  a 

given  frame,  the  term  Ui  has  components 


_ 
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^(Xl(x,<),<) 


dt 


*i(M) 


( 

=  (  fl — 


+  2 


3*ui^.  ,  9ti^  dvk 

5f9xi.  5x1.  dt  5t* 


)■ 


(10.5) 


where  v(x,t)  =  ^Xi(x,t)  =  ^0(t).  The  inertial-type  terms  in  (10.5)  obviously 
occur  because  points  in  Rf  are  moving.  This  also  results  in  the  boundary  conditions 
on  the  boundary  of  Ri'/Si,  where  Si  =  Xi(St,f),  being  with  respect  to  a  moving 
boundary.  Both  of  these  issues  complicate  solving  the  corresponding  boundary  value 
problem,  whether  it  be  by  using  analytical  methods  or  by  constructing  a  finite 
difference  or  a  finite  element  computer  program. 


Fortunately,  the  boundary  value  problem  with  the  balance  of  linear  momentum 
for  phase  2  in  the  form  given  by  Equation  (10.4)  does  not  have  to  be  solved.  Instead, 
Equation  (10.4)  can  be  transformed  into  a  more  tractable  equation.  In  particular,  we 
can  use  the  mapping  x^  and  define  the  function  ui  as 


Ui(x,t)  =  ui(xi(x,  <),<),  (10.6) 

Vx  €  R"*"  and  V<  €  F,  and  then  solve  for  Ui  =  Oi(x,t)  instead  of  Ui  =  Ui(x,,t)  in 
the  boundary  value  problem.  Also,  note  that  Ui(xi,t)  =  Ui(x(xi,t),t)  V  x^  €  Ri' 
at  each  t  €  F. 

In  terms  of  Ui,  the  elastic  potential  and  stress  tensor  for  phase  2  given  by 
Equations  (10.2)i  and  (10.2)2,  respectively,  become 


Wy  =  i  (VuxVtx) .  [Ci(VuiVix)]  -I-  0(|ViU||®)  (10.7) 


ai  =  Ci(Vu»Vtx)-l-0(|Vuin 
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respectively,  Vx  e  R'*'  at  each  *  €  F,  where  Xi  =  Xi(x,t)  in  ViXi. 

Substituting  (10.7)2  into  the  equation  representing  the  balance  of  linear  momentum 
for  phase  2  given  by  (S.9)i  and  expressing  the  acceleration  term  in  that  equation  in 
terms  of  Xi  and  ui  yields,  in  a  given  coordinate  frame. 


dxn  ^Ux^.  dXr^ 
dxfi  dxij  dxtn  3xii 

f  ]  I  f 

^iit‘tydx„dxn  dxi,  dxij  dxm  dxijdxi- j  *' 

J\  dt^  dt^  J' 


(10.8) 


Vx  6  R'*’  at  each  f  €  F,  where  Xi  =  Xi(x,f)  in  V,x.  From  Equation  (10.8),  we 
can  observe  that  not  only  have  most  of  the  inertial-type  terms  been  eliminated,  the 
resulting  boundary  value  problem  is  in  a  completely  Lagrangian  description;  i.e.  it 
is  completely  in  terms  of  the  coordinates  of  the  fixed  reference  configuration.  This 
includes  the  boundary  conditions  for  phase  2  being  specified  with  respect  to  a  fixed 
boundary.  The  only  penalty  that  is  paid  for  this  coordinate  transformation  are  the 
additional  terms  in  Equation  (10.8).  For  the  important  special  case  where  xi  is  a 
homogeneous  deformation,  phase  2  is  a  homogeneous  material  (i.e.  Ci  is  independent 
of  x),  and  fi  =  0,  Equation  (10.8)  given  above  reduces  to 


C. 


dXrn  dXn 
'^^^dxmdxn  dx^  dXi- 


j\  y  ’ 


(10.9) 


where  Xi  =  Xi(x,t)  in  ViX.  For  this  case,  the  additional  terms  in  (10.9)  contribute 
only  to  the  coefficients  of  the  terms  involving  second  spatial  derivatives  of  ui.  Thus, 
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in  at  least  this  case,  the  advantages  of  using  multiple  reference  configurations  along 
with  this  coordinate  transformation  far  outweigh  their  disadvantages. 

For  the  static  problem,  points  in  are  stationary  and  it  is  probably  much 
more  convenient  to  work  in  terms  of  Ui  =  Ui(xi)  instead  of  Ui  =  Ui(x)  in  the 
boundary  value  problem.  In  this  case  the  inhomogeneity  of  Wi  can  be  expressed 
explicitly  in  terms  of  points  Xi  €  R^.  and  the  equation  representing  the  balance  of 
linear  momentum  for  phase  2  would  have  the  same  form  as  the  equation  representing 
the  balance  of  linear  momentum  for  phase  1  given  by  (10.3),  except  that  aU  of 
the  quantities  would  have  1  subscripts  and  the  inertial  terms  would,  of  course,  be 
equal  to  zero. 

For  the  linear  problem  where  it  is  assumed  that  |Vu|  «  1  V  x  € 
R~  and  (Vui|  «  1  Vxi  €  R?"  at  each  f  e  F,  we  can  consider  a  process  where 
the  first  and  second  time  derivatives  of  u  and  ui  are  negligible  Vx  €  R~  and 
V  X  G  R'*',  respectively,  1V|  «  1,  and  is  negligible,  for  all  <  €  F.  Such  a 
process  corresponds  to  a  process  where  the  strains  of  each  phase  are  infinitesimal,  the 
motions  of  each  particle  of  each  phase  relative  to  the  undeformed  configuration  of  that 
phase  are  negligible,  the  magnitude  of  the  phase  boundary  velocity  is  infinitesimal,  and 
the  acceleration  of  the  phase  boundary  is  negligible,  for  all  f  G  F.  In  the  following, 
this  type  of  process  will  be  referred  to  as  a  quasi-static  process,  even  though  it 
does  not  conform  to  the  exact  definition  of  such  a  process.  More  specifically,  this 
type  of  process  is  not  a  true  quasi-static  process  since  time  is  not  just  a  parameter 
in  all  variables  of  the  problem,  and  consequently  the  set  of  all  solutions  as  time 
is  varied  does  not  consist  of  only  static  equilibrium  solutions.  Static  equilibrium 
occurs  only  when  V  =  0.  For  this  quasi-static  process,  it  is  probably  much  more 
convenient  to  solve  for  =  Ui(xi)  instead  of  Ui  =  ui(x)  in  the  boundary  value 


problem.  However,  unlike  the  static  problem,  points  in  Rj"  are  moving.  Therefore, 
the  inhomogeneity  of  Wi  still  needs  to  be  expressed  with  respect  to  points  x  €  R"*" 
(see  Section  6).  Thus,  for  a  quasi-static  process  the  balance  of  linear  momentum  for 
phase  2  would  have  the  same  form  as  it  would  have  for  the  static  case,  except  that 
dCi-.i^j/dxi.  would  be  replaced  by  (dCi..i^^/dx^^(dx„/dxi.y 

11.  The  Linearized  Jump  Conditions 

The  most  direct  form  of  the  conitinuity  of  displacements  condition  is  given  by 
(4.1).  If  we  woik  with  this  condition  in  this  form,  there  is  nothing  to  linearize. 
However,  the  jump  conditions  (4.2)  and  (4.3),  which  together  represent  the  continuity 
of  displacements,  are  such  that  they  or  equations  equivalent  to  them  can  be  linearized. 
In  particular,  as  was  mentioned  in  Section  4,  Equation  (4.2)  is  equivalent  to  Equation 
(4.4).  Since  both  a  and  N  are  real  quantities  in  that  equation,  we  can  conclude  that 
the  two-tensor  a  <g>  N,  and  hence  F;f‘F+  -  F~,  has  two  zero  eigenvalues  and  one 
real  not  necessarily  zero  eigenvalue.^  Thus,  we  can  write  the  characteristic  equation 
for  F^F"*"  -  F“  as 

-A3-h/,A2-/2A  +  /3  =0  (11.1) 


where 

/,  =  tr(F+F+-F-j, 

72  =  i  I  [tr(FrF+  -  F-)]*  -  tr(Fj-F+  -  F-)*|,  (11.2) 

I3  =  det(FtF*  -  F-y 

and 

F+F"*-  -  F-  =  ( Vui)^  +  ( Vu, VixVuO^  -I-  (VuiV.xf  -  ( Vu)".  (1 1.3) 


^  Thii  can  easily  be  seen  by  working  in  the  coordinate  ftame  where  either  N  or  a  coincides  with  a  basis  vector. 
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A  necessary  condition  for  the  existence  of  two  zero  eigenvalues  of  F F'*' — F“  is  that 

/2  =  0,  J3  =  0.  (11.4) 

Equations  (1 1.4)  can  both  be  linearized  with  respea  to  Vu  and  Vui.  The  set  of  all 
vectors  L  that  satisfy  the  linearized  form  of  (4.2)  can  then  be  determined.  Once  this 
is  done,  N(x,  t)  can  be  determined  such  that  N*L  =  0Vx£Sat  each  f  €  F.  N 
can  then  be  substituted  into  the  linear  momentum  jump  condition  given  by  (S.IO) 
yielding  three  scalar  equations  which  can  be  linearized  with  respect  to  Vu  and  Vui. 

(Consider  the  path  of  the  phase  boundary  during  a  motion  in  the  time  interval  F. 
Next  consider  the  special  case  where  we  assume  that  at  each  location  of  the  phase 
boundary  in  this  path  there  is  continuity  of  displacements  at  all  points  on  the  phase 
boundary  for  u  =  0  V  x  e  R~  and  Ut  =  0  V  Xi  6  R*.  In  this  case,  as  can  be  seen 
by  substituting  F~  =  FJ"  =  1  into  (4.4),  F+  must  have  the  form 

F'‘'  =  H-a®n  (11-5) 

at  each  point  on  the  phase  boundary  at  each  t  E  F,  where  the  vectors  a  and  n  may 
both  be  functions  of  x  6  S.^  In  this  case,  we  can  conclude  from  substituting  (11.5) 
into  (4.2)  that  the  first-order  approximation  of  L  is  L,  where  L  •  n  =:  0.  Additionally, 
in  this  case,  the  first-order  approximation  of  the  linear  momentum  jump  condition 
given  by  (5.10)  is 

(ffj  -<7")n  +  p(vJ'-v")K,  =  0,  (11.6) 

where  a  and  <ri  are  given  by  their  respective  linear  constitutive  relations. 

For  the  static  case,  jump  conditions  (4.1)  and  (4.2)  are  the  same,  jump  condition 
(4.3)  is  trivially  satisfied,  and  the  linear  momentum  jump  condition  given  by  (5.10) 

^  This  may  be  a  good  assuffipcion  for  a  material  dut  has  coniimiity  of  displaoemenu  at  the  phase  boaodaty  and 
is  unstressed  when  it  is  in  sutic  equilibrium  with  no  applied  boundary  Iractiont. 
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reduces  to 


(j^JaiF-'^y  -(T-^N  =  0, 


(11.7) 


which  is  equivalent  to  the  continuity  of  traction  across  the  interface.  For  the  case 
where  F'*'  =  1  +  a  (S>  n  (11.7)  reduces  to 

(a+ -<r-)n  =  0.  (11.8) 


For  a  quasi-static  process,  jump  conditions  (4.1)  and  (42)  are  the  same.  However, 
for  this  type  of  process,  the  first-order  approximation  of  jump  condition  (4.3)  is 

If  the  displacements  are  to  be  continuous  at  the  phase  boundary  for  all  t  €  F  with 
respect  to  a  first-order  approximation,  the  mapping  xi  should  be  chosen  such  that 
(11.9)  is  satisfied.  For  the  case  where  the  shape  deformation  is  a  homogeneous 
invariant  plane  strain  for  all  t  €  F,  Xi  must  have  the  form  given  by  (6.9)  in  order  for 
(11.9)  to  be  satisfied  at  each  f  €  F  (see  Section  6).  Also,  for  a  quasi-static  process, 
the  first-order  approximation  of  the  linear  momentum  jump  condition  is  the  same 
as  that  for  the  static  case. 

12.  The  Linearized  Driving  Inaction  and  Kinetic  Reiation 

The  driving  traction  acting  on  the  phase  boundary  was  presented  and  discussed 
in  Section  8.  For  the  linear  case  where  we  assume  that  F'^  =  l  +  a0nat  each  point 
on  the  phase  boundary  at  each  t  €  F,  the  driving  traction  given  by  (8.2)  becomes 
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/  =  f(j  w;y  -  -  i  [(fft  +  a-)a]  ■  a + 0(|va|^  |vb,0  (ixi) 

For  the  quasi-static  case,  using  (4.4)  and  (11.7),  the  driving  traction  given  by  (8.2) 
can  be  written  as 

/  =  (i+W^+  _  W-')  -  a-  •  (f+F+  -  F-).  (12.2) 

For  the  linear  case  and  when  F'*'  =  1  +  a  0  n,  this  driving  traction  becomes 

/  =  {(J -  (^*r)  -  •  a  (12.3) 

where  the  second-order  terms  have  been  neglected.  The  implications  of  this  form  of 
the  driving  traction  for  the  quasi-static  case  will  be  discussed  further  in  Section  13. 

As  mentioned  in  Section  8,  we  can  postulate  a  kinetic  relation  at  the  phase 
boundary  which  relates  the  driving  traction  to  the  normal  component  of  the  phase 
boundary  velocity.  It  is  required  that  this  kinetic  relation  satisfy  (8.4)  so  that  energy 
is  dissipated  (or  conserved  if  the  equality  sign  holds),  instead  of  being  created,  during 
a  martensitic  phase  transformation.  We  next  consider  a  kinetic  relation  which  has  the 
general  form 

>0,  f>h 

=0,  for  fi<f<f2  (12.4) 

.  <0,  /</i 

where  the  constants  /i  and  f2  are  material-dependent  and  ate  such  that  /i  < 


30 


0  and  /2  >  0.^  We  note  that  such  a  kinetic  relation  satisfies  (8.4).  For  the  linear 
case,  a  kinetic  relation  which  has  the  form  given  by  (12.4)  is 


$(/)  = 


0, 


f>f2 

for  /i  <  /  <  /a 
/</i 


(12.5) 


where  f  is  the  linearized  driving  traction  and  vi  and  ui  are  constants  which  depend 
on  the  given  material.  We  note  that  as  i/i  and  U2  -*  oo,  Kt  — f  0  and  the  phase 
boundary  moves  with  the  particles  of  material  at  the  interface  but  does  not  pass  over 
them  converting  them  from  one  phase  to  the  other. 


13.  Reorientation 

The  phenomena  where  a  boundary  traction  is  applied  resulting  in  the  phase 
boundary  separating  two  variants  of  the  same  martensite  moving  and  transforming 
one  variant  into  the  other  is  known  as  reorientation.  In  a  material  that  is  fully 
martensitic,  reorientation  takes  place  until  all  of  the  variants  of  martensite  in  the 
material  are  the  same  variant  of  the  same  martensite,^  or  the  boundary  traction  is 
removed. 

The  issue  of  which  variant  of  martensite  is  piefened  during  the  growth 
process,  whether  the  process  is  reorientation  or  simply  an  austenite-martensite  phase 
transformation  where  several  variants  of  martensite  nucleate  at  different  points  in  the 
material,  is  an  issue  that  has  received  much  attention.  For  the  case  where  a  uniaxial 
tensile  traction  is  ^plied,  it  has  been  observed  from  experiments  in  [12]  that  for 

^  Kinecic  lelitiofu  of  ihi<  fono  have  been  studied  in  [1], 

^  If  this  occurs  in  a  single  crystal  of  austenite,  it  will  transfoim  into  a  single  crystal  of  one  variant  of  marteiisitB 
(see  (12]X 
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18R  martensitic  alloys  the  variant  of  martensite  that  is  preferred  during  the  growth 
process  is  the  variant  which  yields  the  largest  amount  of  extension  due  to  its  shape 
deformation  in  the  direction  of  the  uniaxial  tensile  traction.  The  consistency  of  this 
experimental  observation  with  a  minimum  energy  criterion  is  shown  and  discussed 
in  [4].  There  have  also  been  proposed  criteria  involving  the  shear  traction  on  the 
interface  of  the  martensite,^  and  criteria  based  on  the  shear  stress  on  the  plane  of 
slip  for  internally  slipped  martensites. 

In  the  rest  of  this  section,  the  issue  of  which  variant  of  martensite  is  preferred  for 
the  case  of  a  stress-induced  austenite-martensite  phase  transformation  and  for  the  case 
of  reorientation  for  a  general  state  of  stress  at  the  interface  are  considered.  The  issue 
of  which  variant  is  preferred  for  the  special  case  of  a  uniaxial  tensile  traction  applied 
to  a  cylindrical  body  is  then  investigated.  For  this  problem,  a  result  is  obtained  that 
corresponds  with  observations  made  from  experiments. 

13.1.  The  General  Case 

Consider  the  quasi-static  case  where  the  shape  deformation  of  phase  2  is  a 
homogeneous  invariant  plane  strain.  In  this  case,  the  driving  traction  is  a  special 
case  of  the  driving  traction  given  by  (12.3),  where  a  and  n  are  constant  vectors.  We 
next  assume  that  all  of  the  variants  of  martensite  of  the  material  under  consideration 
have  the  same  value  for  their  elastic  potentials  in  their  undeformed  states,^  and  we 
assume  that  the  material  has  a  kinetic  relation  of  the  form  (12.4). 

To  consider  the  case  of  an  austenite-martensite  phase  transformation,  we  let  phase 
1  represent  the  austenite  and  phase  2  represent  a  variant  of  martensite.  From  the 

”  Such  «  criteria  ii  dtacuued  in.  e.g..  [6], 

"  For  this  comparison  of  the  elastic  potemiala.  they  are  all  considered  to  be  with  respect  to  the  same  reference 
configuratioa 
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assumptions  given  above,  we  can  conclude  that  ^ J  W*j  —  has  the  same 

value  regardless  of  which  variant  of  maitensite  phase  2  represents.  We  next  note  that 
(<T"n)  is  the  traction  on  the  interface  and  that  (c'n)  •  a  is  the  component  of  this 
traction  in  the  direction  of  the  amplitude  vector  of  the  variant  of  maitensite  that  phase 
2  represents  multiplied  by  the  magnitude  of  this  amplitude  vector.  We  also  note  that 
the  variant  that  has  the  largest  negative  value  for  its  driving  traction  will  grow  at 
the  fastest  rate.  From  the  discussion  above,  we  can  conclude  that  this  variant  is  the 
variant  with  the  laigest  value  of  (ff~n)  *  a.  Thus,  for  the  case  under  consideration,  we 
can  conclude  that  the  variant  of  martensite  with  the  lai;gest  value  of  the  component 
of  traction  on  its  interface  in  the  direction  of  its  amplitude  vector  multiplied  by  the 
magnitude  of  its  amplitude  vector  will  be  the  variant  that  is  preferred  during  the 
growth  process  in  a  stress-induced  austenite  to  maitensite  phase  transformation. 

To  consider  the  case  of  reorientation,  we  let  phase  1  and  phase  2  represent  two 
different  variants  of  the  same  martensite.  From  the  assumptions  given  previously, 
we  can  conclude  that  (^J  =  0  for  every  combination  of  variants  that 

phase  1  and  phase  2  can  represent.  Therefore,  the  variant  that  has  an  amplitude  vector 
a  and  an  interface  normal  n  that  result  in  (o'n)  •  a  having  the  largest  value  will 
grow  at  the  fastest  rate.  We  note,  however,  that  in  this  case  a  and  n  for  the  variant 
under  consideration  are  measured  with  respect  to  the  neighboring  variant  across  the 
interface  with  unit  normal  n. 

13,2.  The  Case  of  a  Uniaxial  Tensile  Inaction 

Consider  a  cylindrical  body  parallel  to  the  unit  vector  ei.  Assume  that  a  tensile 
traction  t  =  — <7oei  is  applied  at  the  end  with  unit  normal  — e^,  a  tensile  traction 
t  =  aoei  is  applied  at  the  end  with  unit  normal  ei,  and  the  remaining  surface  of  the 
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cylindrical  body  is  traction  free.  For  this  case,  the  stress  tensor 

r  =  <Toei®ei  (13.1) 

is  a  solution  of  the  field  equations,  the  linear  momentum  jump  condition,  and  the 
boundary  conditions  for  the  quasi-static  process  under  consideration.^  Considering 
(10.1)3,  tl^  first-order  approximation  of  a~  is 

a~  —  (ToCi  ®  ei  (13.2) 

For  this  special  case,  the  linearized  driving  traction  given  by  (12.3)  becomes 

/  =  /er  —  <Toacos^cosa,  (13.3) 

where  far  =  —  {W*)~ ,  a  =  |al,  acos^  =  a  •  ei,  and  cosa  =  n  •  ei. 

We  next  note  that  the  extension  of  a  unit  fibre  of  material  originally  parallel  to 
ei  due  only  to  the  shape  deformation  is 

<5  =  Fei|  -  (eil  =  \J Cj  •  Ce,  -  1,  (13.4) 

where  C  =  F^F.  The  component  of  extension  in  the  direction  of  the  tensile  traction 
is 

6'  =  Cl'  Fci  —  Cl  •  Cl  =  acos^cosa.  (13.5) 

This  is  also  equal  to  the  extension  given  by  (13.4)  if  rotations  are  neglected  or  do 
not  occur.  From  this,  we  can  write  the  driving  traction  given  by  (13.3)  as 

f  =  (13.6) 

Thus,  for  a  given  ao,  the  variant  of  martensite  that  yields  the  largest  component  of 
extension  in  the  direction  of  the  uniaxial  tensile  traction  will  be  preferred  during 

”  This  '•*a  readily  be  seen  by  expreuing  these  equations  in  there  Eulerian  form. 
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the  growth  process.  This  corresponds  with  the  experimental  observations  presented 
in  [12]. 

Let  e  be  a  unit  vector  in  the  (a,  n)  plane  such  that  e  •  n  =  0  and  e  •  a  >  0. 
For  the  case  of  a  uniaxial  tensile  traction  described  above,  the  (shear)  component  of 
traction  on  the  interface  in  the  direction  of  e  is 

5  =  e  •  [(<Toei  0  ei)S]  =  <yocosAcosa  (13.7) 

where  e  •  =  cosA.  For  the  case  where  F  is  a  simple  shear,  e  is  parallel  to  ^ 

which  results  in  A  =  ^,  and  (13.7)  becomes 

S  =  aocos^cosa.  (13.8) 

Thus,  for  this  case,  the  driving  traction  given  by  (13.3)  can  be  written  as 

/  =  /„-aS.  (13.9) 

From  this,  we  can  conclude  that  a  shear  stress  criterion  which  states  that  the  variant 
that  is  preferred  is  the  variant  with  the  largest  value  of  5  will  correspond  to  the 
variant  with  the  largest  negative  driving  traction  for  a  given  ao  only  if  each  variant 
has  a  simple  shear  shape  deformation  and  has  the  same  value  for  |a|. 

14.  An  Applied  Uniform  Hydrostatic  Pressure 

In  this  section,  the  effect  that  an  applied  hydrostatic  pressure  has  on  a  martensitic 
phase  transformation  in  a  given  material  is  studied.^  The  two-phase  elastic  material 
that  was  described  in  the  previous  sections  is  considered.  If  a  given  loading  results 
in  an  increase  in  the  driving  traction,  that  loading  is  considered  to  favor  the  growth 
of  phase  1,  and  if  a  given  loading  results  in  a  decrease  in  the  driving  traction,  that 

”  See  {8]  for  a  discuaaion  of  this  type  of  problem  and  for  a  list  of  some  refetences  where  such  problems  are 
conaideied. 
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loading  is  considered  to  favor  the  growth  of  phase  2.  If,  on  the  other  hand,  the  given 
loading  results  in  no  change  in  the  value  of  the  driving  traction,  neither  phase  is 
favored  by  that  loading.  For  the  problem  under  consideration  it  is  assumed  that  a 
hydrostatic  pressure  exists  such  that 

x  =  -pX  (14.1) 

at  every  point  of  the  deformed  body,  where  p  >  0.  A  quasi-static  process  is 
considered  where  lVul«lVxeR  and  IVujl«lVxj  eR*.  For  this  case,  we  can 

conclude  from  (10.1)3  ^^at  the  first-order  approximation  of  a~  is 

a-  =  -pi  (14.2) 

It  is  assumed  that  the  shape  deformation  of  phase  2  is  homogeneous  and  has  the  form 
F  =  l  +  a<S)n.  For  the  following,  let  (jt/2)  -  8  denote  the  angle  between  a  and  n. 
Note  that  5  is  a  measure  of  the  dilatation  (or  volume  expansion)  of  the  shape 
deformation  of  phase  2.  For  the  case  under  consideration  and  when  5  is  not 
infinitesimal,  the  linearized  driving  traction  given  by  (12.3)  becomes 


/  =  J  w;*  )*  -  (IV*  )■  j  pa  •  n  +  0(l  Vul' ) 
=  f„+  pasin54-0(lVuP,lVjUjl^), 


(14.3) 


Thus,  if  5  is  not  infinitesimal,  we  can  conclude  from  (14.3)  that  a  hydrostatic 
pressure  favors  the  austenite  phase  transformation  in  a  material  where  the  martensite 
phase  is  such  that  sin5  >  0,^'  and  a  hydrostatic  pressure  favors  a  martensite  phase 
transformation  in  a  material  where  the  martensite  phase  is  such  that  sin5  >  0. 
Additionally,  as  can  be  observed  from  (14.3),  these  results  are  independent  of  the 


This  corresponds  with  the  well  known  result  that  a  hydrostatic  pressure  stabilizes  the  austenite 
phase  in  many  iron  alloys  [8].  The  martensite  phases  in  these  iron  alloys  have  small,  but  not 
infinitesimal,  volume  expansions. 
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material  coefficients  of  both  the  austenite  and  martensite.  However,  for  the  case 
where  |<5|  «  1,  the  term  pasin<5  becomes  a  second-order  term.^^  Therefore,  in  this 
case,  the  lowest-order  approximation  of  a  second-order  approximation, 

and  the  second-order  terms  that  have  been  neglected  in  (14.3)  must  be  retained  to 
obtain  a  lowest-order  approximation  of  f  —  /„•  We  note  that  these  second-order 
terms  contain  the  material  coefficients  of  both  phases. 

In  the  rest  of  this  section,  the  case  where  |6|  «  1  is  considered.  The  strains 
in  each  phase  corresponding  to  the  hydrostatic  pressure  are  calculated,  and  the 
continuity  of  displacements  condition  (the  only  nontrivial  portion  of  the  boundary 
value  problem)  is  enforced.  The  driving  traction  is  then  calculated  to  determine 
what  effect  the  hydrostatic  pressure,  the  shape  deformation,  the  material  coefficients, 
and  the  orientation  of  the  phase  boundary  have  on  the  sign  of  f  —  /„'■>  i-e.  on  the 
martensitic  phase  transformation.  The  main  result  of  this  section  is  for  the  case  where 
phase  1  and  phase  2  represent  two  different  variants  of  the  same  manensiie  that  are 
twin  related.  For  this  case,  a  result  that  is  expected  from  physical  considerations  is 
obtained.  These  problems  demonstrate  both  the  convenience  and  accuracy  of  using 
this  continuum  model. 

14.1  The  Assumptions 

Let  {ei,  63,  63}  form  an  orthonormal  basis  for  vectors  in  £^3.  For  the  following 
problem,  a  state  of  plane  strain  is  assumed  with  ea  normal  to  the  plane  of  plane 
strain,  and  a  quasi-static  process  is  considered.  It  is  assumed  that  both  phase  1  and 
phase  2  are  homogeneous  and  have  tetragonal  synunetry.  The  case  of  cubic  symmetry 
for  phase  1  and/or  phase  2  can  be  considered  as  a  special  case,  once  the  general 

’’  Note  that  because  v  in  (14.3)  is  related  to  the  infinitesinial  strains  through  the  constitutive  equations,  it  is 
consideied  to  be  a  first-order  term  in  the  equation  representing  the  driving  traction. 
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results  are  obtained.  It  is  assumed  that  each  phase  has  preferred  directions  [13]. 
For  the  case  considered  here  where  the  tetragonal  symmetry  reflects  the  synunetry 
of  a  tetragonal  crystal  lattice,  the  three  preferred  directions  of  each  phase  can  be 
considered  to  be  in  the  [1  0  0],  [0  1 0],  [0  0  1]  directions  of  a  tetragonal  unit  cell  of 
each  phase,  respectively.  Let  {hj",  h+,  h,  }  represent  the  preferred  directions  of  phase 
2  and  assume  that  {ej,  ej,  63}  are  such  that  they  are  aligned  with  {hj‘,hj,hj}, 
respectively.  Let  hj,  hj}  represent  the  preferred  directions  of  phase  1,  with  hj 
in  the  direction  of  63  and  the  angle  between  and  ei  denoted  by  a(Figure  3).  For 
this  plane  strain  problem,  it  is  assumed  th^  the  shape  deformation  is  homogeneous 
with  F  =  l+a(Sn,  where  a  and  n  are  in  the  plane  of  plane  strain;  i.e.,  in  the  {e^,  63} 
plane.^^  Let  <f>  denote  the  angle  between  n  and  ei,  and  let  S  be  defined  as  it  was  in 
the  beginning  of  this  section.  As  mentioned  above,  it  is  assumed  that  |6|  «  1  in  the 
following.  It  is  also  assumed  that  a  plane  strain  hydrostatic  pressure  exists  such  that 

T  =  —pi,  (14.4) 


at  every  point  of  the  deformed  body.  Additionally,  in  the  following  analysis,  we 
will  decompose  the  displacement  gradients  into  their  synunetric  and  skew-syirunetric 
parts  as  follows 


Vu  =  c  +  w, 


ViUi  =€1 


(14.5) 


where  (u  =  skewVu,  wi  =  skewViUi. 


Unleu  oifaerwue  indicated,  for  the  following  plane  itiain  problem,  it  is  assumed  that  all  tensors,  including 
vecttxB,  correspond  to  the  two-dimensional  Euclidean  space  containing  (ci,  ot). 
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14^  The  Stress  and  Infinitesimal  Strain  Tensors 


Considering  (14.4),  (10.1)3,  and  (10.2)3,  we  can  conclude  that  the  first-order 
approximations  of  a  and  ai  are 


a  =  -pi. 


(Ti  =  -pi, 


(14.6) 


respectively.^ 

For  phase  2,  because  {hjf',  h^}  coincide  with  {ei,  63),  the  inplane  components 
of  infinitesimal  strain  tensor  Ci  in  the  {e^,  63}  frame  as  functions  of  the  components 
of  stress  in  that  frame  are 


'Ktni 

^U22 

0 

1  £+ 

»  = 

0 

I  e+ 

V  ^1,  J 

0 

0 

(14.7) 


where  denote  the  components  of  the  compliance  four-tensor  for  phase  2  in  the 
{ei,  63}  frame  (see  [3],  [7]).^®  Substituting  (14.6)2  into  (14.7)  yields 


'  Ktm  +  1 

»  =  -p< 

‘  > 

0  J 

(14.8) 


For  phase  1,  the  inplane  components  of  infinitesimal  strain  e  in  the  {e^,  63} 
frame  in  terms  of  the  stress  given  by  (14.6)i  are 


^  Note  that  tinoe  r.  ei.  t,  and  «t  are  homogeneoua,  =  v,  et  =  ri,  e~  =  a,  and  =  «i. 

Note  that  the  compliance  four-tensor  is  the  inverse  of  the  elasticity  four-tensor.  Also,  in  [3]  and  fT].  oithotropic 
symmetry  is  considefed.  Orthotropicsymmetry  and  the  type  of  tetragonal  symmetry  considered  here  ate  equival^ 
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I{^riii(cosaf  +  A:^22(sina)^  +  K^x-a)  1 

{■^llllC*^®)  ■^'  -^1122}  I  ’ 

K-^1111  ~  K.^)siii2a  j 

where  K~0^g  denote  the  components  of  the  compliance  four-tensor  for  phase  1  in 
the  {hj",  hj}  frame. 

143.  The  Crntinuity  of  Displacement  Condition 

Because  the  true  stress  is  assumed  to  be  uniform  and  constant  throughout  the 
body,  the  linear  momentum  field  equations  and  the  traction  jump  condition  are 
trivially  satisfied.  A  necessary  condition  for  the  continuity  of  displacements  at  the 
phase  boundary  is  that 

det(F+F+-F-)  =0.  (14.10) 

For  the  linear  case  that  is  considered  here,  the  first-order  ^proximation  of  L  is  L, 
where  L  •  n  =  0.  In  terms  of  the  displacement  gradients,  the  above  equation  becomes 

det(a®  n-H  Viuj'(a0n)  +  —  Vu”)  =  0.  (14.11) 

Using  (14.5),  the  first-order  approximation  of  the  above  equation  can  be  written  as 

-  eri)sin2<^  +  er,cos2<^  +  w”  -  =  0 ,  (14.12) 

where  lj~  and  u*  denote  the  and  components  of  u  and  in  the  {ei,  63} 
frame,  respectively,  and  has  been  taken  to  be  equal  to  zero  because  of  (14.8). 
We  note  that  for  the  linear  case  under  consideration,  the  symmetric  parts  of 
the  displacement  gradients  of  u  and  ui  are  completely  determined  by  the  stress 
distribution  and  are  given  by  (14.8)  and  (14.9),  respectively,  and  Equation  (14.12) 
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is  the  only  equation  that  restricts  the  skew-symmetric  parts  of  these  displacement 
gradients.  Substituting  (14.8)  and  (14.9)  into  (14.12)  and  solving  foru;'^  — a;~  yields 

u;"*"  —  u>  = —(Jf^22  ~  +  2^-^2222  “  —  2^) .  (14,13) 


14.4.  The  Driving  Traction 

The  second-order  approximation  of  the  driving  traction  corresponding  to  the 
problem  under  consideration  is  as  follows: 

f  ~  fcr+p{a6 +  h>  [(ViUi)a)  -|-  i trei  -  i trc}  ,  (14.14) 

where  a  and  h  are  given  by  their  lowest-order  approximations.  Using  (14.5),  (14.8), 
and  (14.9),  the  second-order  approximation  of  the  driving  traction  given  above  can 
be  written  as 

/  =  /er  +  p{a5  +  auj'*'  +  i  ap(iiL’i'^ni  “  K^)sw2<^  -I-  5p[tr(Kl)  -  tr(Kil)]}  . 

(14.15) 


14.5.  The  Nondimensional  Form  of  the  Driving  Traction  and  Infinitesimal 
Rotations 

Equations  (14.13)  and  (14.15)  can  be  written  in  nondimensional  form  by  using 
the  following  nondimensional  quantities 


2au;" 

“  pxtm ' 

. .  2au;'*' 

~PKW 

II 

-  _  Kfiii 

-^2222 

^2  ~ 

II 

2/ 

1  j^+  » 

(14.16) 
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/cr  = 


2/c 


6  = 


2a6 


nil 


pKtx 


nil 


Kt. 


[tr(Kl)-(K,l)]. 


nil 


The  nondimensional  form  of  Equation  (14.13)  is 

—  l^sin2^  +  —  ^f^sin(2a  —  2<i>) .  (14.17) 

The  nondimensional  driving  traction  /  is 

/  =  /cr  +  ^+w++4  +  a(l-X^)siii2^.  (14.18) 

Note  that  is  a  parameter  in  the  above  equation  for  /.  We  can  obtain  /  as  a 
function  of  Cj"  by  using  (14.17).  The  resulting  equation  is 

/ss  <^(^2  •“  ^f)8ui(2a  —  2^) .  (14.19) 


14.6.  An  Austenite-Martensite  Phase  Dransformation 

To  investigate  the  effect  of  the  hydrostatic  pressure  on  an  austenite-maitensite 
phase  transformation,  we  assume  that  phase  1  represents  the  austenite  and  phase 
2  represents  the  martensite.  From  the  discussion  in  the  beginning  of  this  section, 
a  hydrostatic  pressure  will  favor  the  austenite  phase  transformation  if  it  results  in 
/—  /er  >  0,  and  a  hydrostatic  pressure  will  favor  the  martensite  phase  transformation 
if  it  results  in  /—/«.<  0.  We  note  that  a  value  for  ui"  and  in  (14.18)  and 
(14.19),  respectively,  might  be  prescribed  by  some  displacement  boundary  condition 
that  is  applied  to  the  material  in  addition  to  the  hydrostatic  pressure  and  that  is 
consistent  with  the  assumed  state  of  stress  within  the  material.  Thus,  from  (14.18), 
(14.19),  and  the  nondimensional  variables  given  by  (14.16),  one  can  observe  what 
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values  for  the  shape  deformation  of  the  maitensite,  the  material  coefficients  of  both 
phases,  the  orientation  of  the  phase  boundary,  and  the  additional  boundary  conditions 
will  result  in  the  hydrostatic  pressure  favoring  the  martensite  phase  transformation 
and  what  values  of  these  quantities  will  result  in  the  hydrostatic  pressure  favoring 
the  austenite  phase  transformation.  Also,  we  note  that  as  6  increases  in  magnitude, 
the  results  that  were  obtained  above  for  the  second-order  case  {q)proach  the  results 
that  were  obtained  in  the  beginning  of  this  section  for  the  first-order  case,  as,  of 
course,  should  be  expected. 

14.7.  The  Case  of  Twin  Related  Variants 

Consider  the  case  where  phase  1  and  phase  2  represent  two  variants  of  the  same 
maitensite  which  are  twin  related  and  the  isothermal  process  is  at  a  temperature  near 
the  transformation  temperature.  In  this  case  we  have:  (1)  £  =  0,  since  one  twin 
has  a  simple  shear  shape  deformation  relative  to  the  other  twin;  (2)  fcr  =  0,  since 
JWi{l)  =  IV(1)  for  two  variants  of  the  same  martensite  near  the  transformation 
temperature;  (3)  since  both  phases  represent  the  same  material;  (4) 

^  =  0,  as  a  result  of  (3);  and  (S)  a  =  2<^,  as  can  be  seen  from  Figure  3.2.  In  this 
case,  the  driving  traction  given  by  (14.18)  becomes 

/  =  -I-  a  (l  -  sin2<i> ,  (14.20) 

and  Equation  (14.17)  becomes 

-  w-  =  -2a(l  -  k^yui2<f>.  (14.21) 

If  we  now  require  that  the  deformation  of  one  twin  be  symmetric  with  respect  to 
the  deformation  of  the  other  twin  (due  to  the  symmetric  stress  distribution),  we  must 
require  that 


w”  =  — 


(14.22) 
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In  this  case,  substituting  (14.22)  into  (14.21),  solving  for  in  the  resulting  equation, 
and  then  substituting  the  resulting  equation  for  u'*'  into  (14.20),  yields 

/  =  0.  (14.23) 

This  is  what  should  be  expected,  since  everything  else  in  the  problem  is  symmetric, 
and  consequently  one  twin  should  not  be  preferred  over  the  other. 

15.  Summary  and  Concluding  Remarks 

In  this  paper,  a  continuum  model  for  materials  that  can  undergo  martensitic  phase 
transformations  was  developed.  The  continuum  model  was  then  used  to  study  several 
problems  that  deal  with  which  phase  or  which  variant  of  martensite  is  preferred 
during  the  application  of  a  mechanical  loading. 

In  the  continuum  model  that  was  developed  in  this  paper,  each  phase  has  its  own 
constitutive  equations  which  are  defined  with  respect  to  a  reference  configuration  that 
coincides  with  an  unstressed  undeformed  configuration  of  that  phase.  This  reference 
configuration  is  also  used  for  the  expression  of  the  field  equations  for  that  phase.  Widi 
respect  to  the  undeformed  austenite  phase,  these  reference  configurations  coincide 
the  shape  deformations  of  the  martensites.  The  field  equations  and  juii^>  conditions 
for  a  two-phase  material  in  a  purely  mechanical  process  were  derived  and  discussed. 
The  general  form  of  the  constitutive  equations  that  are  most  appropriate  for  the  case 
where  each  phase  of  the  material  behaves  elastically  were  presented  and  discussed. 
It  was  pointed  out  that  the  mapping  representing  the  reference  configuration  for  each 
phase  should  be  such  that  the  stress-power  for  that  phase  is  equal  to  the  time  rate  of 
change  of  the  integral  of  the  elastic  potential  for  that  phase,  where  the  elastic  potential 
and  stress  tensor  for  that  phase  are  defined  with  respect  to  the  reference  configuration 
for  that  phase.  This  is  the  case  if  the  mapping  for  each  phase  has  the  form  given 


by  (6.7).  Additionally,  one  can  argue  that  for  a  purely  mechanical  process  and  a 
material  with  phases  that  behave  elastically  xi  should  always  have  the  form  given 
by  (6.7)  since  this  is  the  form  that  xx  must  have  in  order  for  F,  and  consequently 
the  transformation  strain,  to  be  independent  of  time.  It  was  also  pointed  out  that 
if  the  reference  configuration  for  a  given  phase  is  moving,  the  quantities  that  are 
defined  with  respect  to  that  reference  configuration  (e.g.  the  stress)  do  not  have  a  true 
nominal  form.  Because  different  elastic  potentials  are  used  for  the  different  phases 
of  the  material,  the  material  symmetry  groups  of  the  different  phases  can  be  chosen 
independently  of  each  other  and  can  be  chosen  to  reflect  any  type  of  crystal  symmetry 
with  any  orientation.  Thus,  the  change  in  crystal  structure  that  takes  place  during 
a  martensitic  phase  transformation  and  the  orientation  relation  between  the  crystal 
lattices  of  the  different  phases  of  the  material  can  be  accurately  represented.  In  fact, 
because  the  elastic  potential  for  each  phase  is  defined  with  respect  to  the  undeformed 
configuration  of  that  phase,  the  material  syiiunetry  group  for  that  elastic  potential 
can  be  chosen  to  be  the  crystallographic  point  group  corresponding  to  the  crystal 
sjrmmetry  of  that  phase.  The  driving  traction  and  kinetic  relation  corresponding  to 
the  two-phase  material  were  also  considered. 

One  of  the  main  advantages  of  working  with  this  continuum  model  is  that  the 
field  equations  and  jump  conditions  are  in  forms  that  permit  direct  linearization, 
while  still  retaining  finite  shape  deformations  for  the  martensites.  This  is  the  case 
since  the  field  quantities  for  each  phase  are  in  terms  of  displacements  that  ate 
measured  from  the  reference  configuration  corresponding  to  an  unstressed  undeformed 
configuration  of  that  phase,  and  consequently,  for  the  appropriate  initial  and  boundary 
conditions,  the  displacement  gradients  can  be  considered  infinitesimal.  The  linearized 
field  equations,  jump  conditions,  driving  traction,  and  kinetic  relation  were  then 
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considered.  It  was  pointed  out  that  working  with  displacements  that  are  functions  of 
points  on  a  moving  reference  configuration  results  in  the  accelerations  in  terms  of 
these  displacements  containing  many  inertial-type  terms  and  results  in  the  boundary 
conditions  for  the  phases  having  to  be  expressed  with  respect  to  moving  boundaries.  A 
coordinate  transformation  that  eliminates  most  of  these  inertial-type  terms  and  results 
in  a  completely  Lagrangian  description  of  the  boundary  value  problem  while  still 
wotidng  with  field  quantities  that  are  defined  with  respect  to  the  different  reference 
configurations  was  discussed.  Another  advantage  of  using  this  continuum  model  is 
that  for  the  linear  problem  the  true  stress  for  each  phase  is  approximately  equal  to 
the  nominal  stress  for  that  phase.  This  is  very  useful  for  solving  certain  types  of 
boundary  value  problems,  including  the  linear  problems  that  were  considered  in  the 
rest  of  this  paper.  All  of  these  problems  dealt  with  the  issue  of  which  phase  or  which 
variant  of  martensite  is  preferred  during  the  growth  process  during  the  application  of 
a  mechanical  loading.  After  the  general  problem  was  discussed,  the  special  case  of 
an  applied  uniaxial  tensile  traction  was  considered.  A  result  that  corresponds  with 
observations  that  were  made  from  experiments  was  derived.  The  last  problem  that 
was  considered  in  this  paper  involves  the  application  of  a  hydrostatic  pressure  to  a 
two-phase  material  with  a  martensite  phase  that  has  a  finite  shape  deformation  with  an 
infinitesimal  dilatation.  The  case  where  the  two  phases  represent  two  different  variants 
of  the  same  martensite  was  then  considered  as  a  special  case.  It  was  shown  that  ftn* 
this  case  the  hydrostatic  pressure  favors  the  growth  of  neither  variant,  whidt  is  what 
should  be  expected  since  everything  else  in  the  problem  is  completely  symmetric. 
All  of  the  problems  that  were  considered  in  this  paper  demonstrate  the  convenience 
and  accuracy  of  using  the  continuum  model  that  was  developed  in  this  paper. 
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